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In this letter, we revise the QED gauge invariance for the hadron tensor of Drell-Yan type processes with 
the transversely polarized hadron. We perform our analysis within the Feynman gauge for gluons and make a 
comparison with the results obtained within the light-cone gauge. We demonstrate that QED gauge invariance 
leads, first, to the need of a non-standard diagram and, second, to the absence of gluon poles in the correlators 
{ilyy±A + ip) related traditionally to dT(x, x) / dx. As a result, these terms disappear from the final QED gauge 
invariant hadron tensor. We also verify the absence of such poles by analysing the corresponding light-cone 
Dirac algebra. 
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I. INTRODUCTION 


II. KINEMATICS 


In the recent times, we have observed the renaissance in 
the nucleon structure studies through the Drell-Yan type pro¬ 
cesses in the existing (FermiLab, Relativistic Heavy Ion Col¬ 
lider, see Hl|2l) and future (J-Parc, NIC A) experiments. One 
of the most interesting subjects of such experimental studies 
in this direction is the so-called single spin asymmetry (SSA) 
which is expressed with the help of the hadron tensor, see for 
instance El or 00 . 

Lately, we have reconsidered 0 this process in the contour 
gauge. We have found that there is a contribution from the 
non-standard diagram which produces the imaginary phase 
required to have the SSA. This additional contribution leads 
to an extra factor of 2 for the asymmetry. This conclusion 
was supported by analysis of the QED gauge invariance of the 
hadron tensor. 

In comparison, the analysis presented in (7) which uses the 
axial and Feynman gauges does not support the latter conclu¬ 
sion. For this reason, we perform here the detailed analysis of 
hadron tensor in the Feynman gauge with the particular em¬ 
phasis on the QED gauge invariance. We find that the QED 
gauge invariance can be maintained only by taking into ac¬ 
count the non-standard diagram. Moreover, the results in the 
Feynman and contour gauges coincide if the gluon poles in 
the correlators ( , ip'y±A + 'il>} are absent. This is in agreement 
with the relation between gluon poles and the Sivers function 
which corresponds to the "leading twist" Dirac matrix 7 +. We 
confirm this important property by comparing the light-cone 
dynamics for different correlators. 

As a result, we derive the QED gauge invariant hadron ten¬ 
sor which completely coincides with the expression obtained 
within the light-cone contour gauge for gluons, see 0 . 


We study the hadron tensor which contributes to the single 
spin (left-right) asymmetry measured in the Drell-Yan process 
with the transversely polarized nucleon (see Fig. [I]): 

N^\ Pl ) + N{p 2 ) -4 7 *{q) + X(P x ) 

t{h)+t{h) + X{P x ). ( 1 ) 


Here, the virtual photon producing the lepton pair (h+h = (?) 
has a large mass squared ( q 2 = Q 2 ) while the transverse mo¬ 
menta are small and integrated out. The left-right asymmetry 
means that the transverse momenta of the leptons are corre¬ 
lated with the direction Sxe, where ,5), implies the transverse 
polarization vector of the nucleon while e z is a beam direction 
0 . 

Since we perform our calculations within a collinear fac¬ 
torization, it is convenient to fix the dominant light-cone di¬ 
rections as 


Q * Q 

Pl * x B V2 H ’ ?2 "W2 
n*" = (1/V2, 0 T , 1A/2), n 


n, 

M _ 


( 2 ) 

(1/V2, Or, —l/v/2). 


So, the hadron momenta pi and p 2 have the plus and minus 
dominant light-cone components, respectively. Accordingly, 
the quark and gluon momenta k\ and i lie along the plus di¬ 
rection while the antiquark momentum k 2 - along the minus 
direction. The photon momentum reads (see Fig. [TJ 

<7 = l\ + I 2 = k\ + k 2 (3) 


which, after factorization, will take the form: 


q = Xipi + yp 2 + q T - 


(4) 


III. THE DY HADRON TENSOR 
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Fig. [TJ the left panel) reads 

dWfs* nd .) = J d 4 k i d 4 fc 2 <5 (4) (fci +k 2 -q)x 

J d A e<f> i a )bl3] (k 1 ,e)$ [ ' y - ] (k2) x 

tr[7 /i 7 /3 7 y 7 + 7 Q S'(f - k 2 )\ , 


where 


4 A) 17/51 (k u l) = 


d ~2 


(pi, s T \i’{mh 0 gA a (z)ip{o)\s T ,p ± ) 


< £ >[7 ] (k 2 )=Fi (p 2 |V’(t?2)7 ^(0)|P2 


(5) 


( 6 ) 


(7) 


Throughout this paper, T\ and T 2 denote the Fourier transfor¬ 
mation with the measures 


d 4 y 2 e ifc2 '" 2 and d 4 r h d 4 ze~ ikl ' ril - u ' z , 


( 8 ) 


respectively, while T x 1 and T 2 1 mark the inverse Fourier 
transformation with the measures 


dy e iyX and d Xl dx 2 e ixiXl+i ^~ x ^ x \ 


(9) 


We now implement the factorization procedure (see for in¬ 
stance mm ) which contains the following steps: (a) the 
decomposition of loop integration momenta around the cor¬ 
responding dominant direction: ki = Xip + (ki ■ p)n + kx 
within the certain light cone basis formed by the vectors p 
and n (in our case, n* and ?r); (b) the replacement: d 4 ki => 
d 4 ki dxiS(xi — ki ■ n) that introduces the fractions with the 
appropriated spectral properties; (c) the decomposition of the 
corresponding propagator products around the dominant di¬ 
rection. In Eqn. <0, we have (here, Xij = Xi — Xj) 


S{1 - k 2 ) = S(x 2 ipi - yp 2 ) + 


( 10 ) 


dS(£ - k 2 ) 

d£o 


k 2 —yp 2 


t T 


e=x 2 ipi 

(d) the use of the collinear Ward identity: 

dS(k ) 


dk n 


= S(k)'y p S(k), S(k) = 


-t 


k 2 


- IE 


(e) performing of the Fierz decomposition for ip a (z) 4>p(0) in 
the corresponding space up to the needed projections. 

After factorization, the standard tensor, see Eqn. < 0 , is split 
into two terms: the first term includes the correlator without 
the transverse derivative, while the second term contains the 


correlator with the transverse derivative, see Eqns. (IOi and 

©-©■ 

The non-standard contribution comes from the diagram de¬ 
picted in Fig. [T] (the right panel). The corresponding hadron 
tensor takes the form 0 : 


dVVftj' , i \ 

(Non-stand.) 


(ii) 


d 4 k x d 4 k 2 c) (4 )(fci + k 2 - q)tT[y p J r (ki)Y^{k 2 )\ , 


where the function J-{k\) reads 

T(k x ) = 5(fcr)7“ J d 4 r/i e ~ iki m x 

(j> x , S T \f>(r ll ) gA a (0) ip(0)\S T ,p x ) . 


( 12 ) 


For convenience, we introduce the unintegrated tensor 
for the factorized hadron tensor W Ml/ of the process. It reads 

W LV = J d 2 q T dW pu = ^ J d 2 q T S^ 2 \q T ) x 
i J dxi dy [5 (xi/x b - l)S(y/y B - (13) 

After calculation of all relevant traces in the factorized hadron 
tensor and after some algebra, we arrive at the following con¬ 
tributions for the unintegrated hadron tensor (which involves 
all relevant contributions except the mirror ones): the standard 
diagram depicted in Fig. [T] the left panel, gives us 


W 


(Stand.) 


Pl e vS T -p2 


W (Stand., 9 l) — Q{y) | (14) 

dx 2 -—- °^—B ( ' 1 \x 1 ,x 2 ) 

B { - 2 \x 1 ,X 2 ) 


y J x\ - x 2 +ie 

x\ I dx 2 


P2 £ pS T - P 2 + P2 L£ vS T -p2 

LXi X\ 


Xi - x 2 + ie 


Pl £ uS T — P 2 


d X ‘ 


B ( ' JL \x !,X 2 ) 
Xi - x 2 + ie 


while the non-standard diagram presented in Fig. |TJ the right 
panel, contributes as 


W 


(Non-stand.) did) ^ 
X\ 


P* 2 _ vS T —p2 


J dx 2 {B (l \ Xl ,x 2 ) + B ( 2 \xi,x 2 )y (15) 

Here we introduce the shorthand notation: £ abcd _ 
P 2 P 3 P 4 A tJ i B ll2 C fl3 D ll4 with e 0123 = 1. Moreover, the 
parametrizing functions are associated with the following cor¬ 
relators: 


is a+sT (pip 2 ) B^\xi,x 2 ) = 


B 2 


(jpi, s T \ip(vi) 7 + gA°i{z) v>(o)|s ,T ,pi) 


i£ +l3sT {p l p 2 )B ( - 2 \xi,x 2 ) = 


d ~2 


{pi, S T \'tp(r)i) 7 ^ gA + (z) ip(0)\S T ,p 1 ) 


ipfe p+sT (pip 2 )BA)(xi,x 2 ) = 

(pi,S T \f)(r) 1 )i + ( d p ± gA + (z )) •0(O)|5' T ,pi) 


(16) 


(17) 


(18) 


where rji = Ain, z = A 2 n, and the light-cone vector ft is a 
dimensionful analog of n ( n~ = pf /{piP 2 ))- 
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As known from J 6 ), the function B ,yV) [x \. x 2 ) for the DY 
process can be unambiguously written as 


B < ' 1 \x 1 ,x 2 ) 


T{x !,X 2 ) 

Xi — X 2 + ie ’ 


(19) 


where the function T(x i,x 2 ) £ Ke parametrizes the corre¬ 
sponding projection of (ip G a p ip), i.e. 


Notice that we have derived (see j 6 |) the certain complex pre¬ 
scription in the r.li.s. of ( [T9| within the contour gauge. In 
this letter, we assume that the same prescription takes place 
in the Feynman gauge too lfl3l . With respect to the functions 
(xi, x 2 ) and B^ ± \x i,x 2 ), we demonstrate below that 
these functions do not possess the gluon poles and, therefore, 
cannot be presented in the form of (fl9|). 


£ a+ sT (p lP2 )T(Xi,X 2 ) = 


T 2 


(pi,S T \ip{m) 7 + n u G%*(z) ip(0)\S T ,pi) 


( 20 ) 


Summing up all contributions from the standard and non¬ 
standard diagrams, we finally obtain the expression for the 
unintegrated hadron tensor as 


^ = WfsU) + Wfstand., + 1 ^(Non-stand.) = «(») 


-AAV 


-AAV 


-AAV 


Pi £ vS T — P2 


U p „ 

P%_Pi_ 
lzi y - 


P^e vST [dx 2 BW(x 1 ,x 2 )- rP 2 £ MS T - P2 + P f L £ vS T - P2 '\x 1 [ 

Xi J IXi Xi J J 

B <y± -\x 1 ,x 2 ) 1 


£ vS T -p 2 J dx2 Bi 1 \ Xl ,x 2 ) + 
-B ( 2 ) (xi,x 2 ) , 


dx 


J Xi - x 2 + ie 


dx 2 


‘ xi — x 2 + ie 


( 21 ) 


Notice that the first term in Eqn. 0 coincides with the 
hadron tensor calculated within the light-cone gauge A + = 0. 


IV. QED GAUGE INVARIANCE OF HADRON TENSOR 


Let us now discuss the QED gauge invariance of the hadron 


tensor. From Eqn. (21 1 , we can see that the QED gauge in¬ 
variant combination is 


_ 


_ LA IA 

P2__Pl_ pV S T - P2 
Lxi y 1 

with = ?„7 7 "'= 0. 


( 22 ) 


We can see that there is a single term with p 2 which does not 
have a counterpart to construct the gauge-invariant combina¬ 
tion 


la 

P2_ 

X\ 


LA 

P __[ 

y ' 


(23) 


Therefore, the second term in Eqn. ( fl4| ) should be equal to 
zero. This also leads to nullification of the second term in 


Eqn. ([13 1 . 

Hence, the only way to get the QED gauge invariant com¬ 
bination (see ( |22| ) is to combine the first terms in Eqns. ( p~4] > 
and CD- This combination justifies the treatment of gluon 
pole in x 2 ) using the complex prescription. 

In addition, we conclude that the third term in ( [T4| does not 
contribute to SSA. 

The suggested proof explores only the gauge and Lorentz 
invariance. Let us consider the other reasoning to justify these 


properties of correlators, starting with the correlator which 
generates the function B^ 2 ' (xi,x 2 ): 

J (, d\ 1 dX 2 )e - ixiXl - i ^- Xl '> x2 x 

(Pi, S T \ip(X 1 h) 7 ^ A+(A 2 n) ^(0)|S T ,pi) = 

ie +l3sT - {pip 2 ) B (2) (x i,x 2 ) ■ (24) 

We are going over to the momentum representation for the 
correlator from the l.h.s. of Eqn. ([24|. Schematically, we have 


u(k 1 )'Ypu{k 2 ) 


1 


X .... X 


P + ie ’ 


(25) 


where the gluon momentum is £ = k 2 — k\ and k\ = 
(x\Pi , kp), kuj, k 2 = (x 2 pi,k 2 ,k 2 ±). This situation has 
been illustrated in Fig. [2] see the left panel. Up to the order 
of g, we are also able to write down that (see Fig. [2| the right 
panel) 


u(k 1 )'yjS(k 2 )u(k 1 ) 


1 


x .... x 


P + ie 


(26) 


where S(k 2 ) = S(k 2 )"f + . From both these equations, it is 
clear that to get the non-zero contribution we must have either 
ki j_ 7^ 0 or k 2 _l 7^ 0. Indeed, 


(fci)7^5(fc 2 )u(fei) => Spki+k! = k 2f) kj + kfpkt (27) 


Therefore, the gluon propagator in Eqns. (25 1 and (26 1 takes 
the following form (cf. flOl '): 


1 


P + ie 2(x 2 — xi)pf(~ — f \+ie 


( 28 ) 
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One can conclude that, in the case of the substantial trans¬ 
verse component of the momentum, there are no sources 
for the gluon poles at x \ = x 2 . As a result, the function 
B^ 2 \x i , x 2 ) has no gluon poles and, due to T-invariance CD 
(B ( '' 2 \x\,X 2 ) = —B( 2 \x 2 ,xi)), obeys B^ 2 \x,x) = 0. 

On the other hand, if we have 7 + in the correlator (see Eqn. 
©), the transverse components of gluon momentum are not 
substantial and can be neglected. That ensures the existence 
of the gluon poles for the function B^\xi, X 2 ). This cor¬ 
responds to the fact that the Sivers function, being related to 
gluon poles, contains the "leading twist" projector 7 +. More¬ 
over, we may conclude that the structure 7 + (<T L A + ) does not 
produce the imaginary part as well as SSA in the Feynman 
gauge. 


V. CONCLUSIONS AND DISCUSSIONS 


Working within the Feynman gauge, we derive the QED 
gauge invariant (unintegrated) hadron tensor for the polarized 
DY process: 


r-r-rUl/ t-tt/TI/ 

Wqi = W( Non _ stand ) + W( Stand ) 


q(y) 


P2_ 

LXi 


U 

Pi 


£ vS p 2 dx 2 B^ (x 1 ,x 2 ). (29) 


After calculating the imaginary part (or, in other words, after 
adding the mirror contributions), and, then integrating over x d 
and y (see Eqn. (| 1 3}), we get the QED gauge invariant hadron 
tensor as 


Wg = q(y B ) 


u 

Xb 


U 

Pl 

VB 


~vS ~P2 


T(x b , x B ) ■ (30) 


This expression fully coincides with the hadron tensor which 
has been derived within the light-cone gauge for gluons. 

Moreover, the factor of 2 in the hadron tensor that we found 
within the axial-type gauge m is still present in the frame of 
the Feynman gauge. In order to show this factor of 2, let us 
introduce the mutually orthogonal basis (see ED) as 


Zp = Plp -P2p = X B Plp ~VBP2p (31) 

and 

^ p = — ^ PP 1 P 2 Q • (32) 

Here p) p are the partonic momenta (q^ = pi p + P 2 p)- With 
the help of ( |Xi~[ ) and ( [32] ), the lepton momenta can be written 
as (this is the lepton c.m. system) 

hp=\qp + jfp(e,ip-,X,Y,Z), 

l2p=\qp-®fp(e,y,X,Y,Z), (33) 

where A = A/\/—A 2 and 

fp(9,(p-,X,Y,Z) = (34) 

Xp cos ip sin 8 + Yp sin <p sin 9 + Zp cos 9 . 



Within this frame, the contraction of the lepton tensor with the 
gauge invariant hadron tensor ( [30] ) reads 

Cp V Wg = -2cos 6£ hsT P^q(y B )T{xB,XB). (35) 

We want to emphasize that this expression in ( [35} differs by 
the factor of 2 in comparison with the case where only one 
diagram (presented in Fig. [I] the left panel) has been included 
in the (gauge non-invariant) hadron tensor, i.e. 

^W^ ±) = \Cp V W^. (36) 

Therefore, from the practical point of view, if we neglect the 
diagram in Fig. |T] (right panel) or, in other words, if we use 
the QED gauge non-invariant hadron tensor, it yields the error 
of the factor of two. 

Further, based on the light-cone dynamics we argue that 
there are no gluon poles in the correlators (’ij)'y±A + 'ip). This 
means that the function B^(x 1 , x 2 ) does not have the repre¬ 
sentation similar to m We also show that the Lorentz and 
QED gauge invariances of the hadron tensor calculated within 
the Feynman gauge require that the function B (2> (xi, x 2 ) can¬ 
not have gluon poles. 

The fact that the function B ( ' 2> (xi, x 2 ) cannot be presented 
in the form of ( | 1 9} directly leads to the absence of dT/dx 
in the final expression of the gauge-invariant hadron tensor. 
Indeed, from (|14j), one can see that B {1> (x -\, . 7 : 2 ) contributes 
to the standard hadron tensor as 


p ^ ST -P2 + p^S T ~P2 


5 (2) (x i,x 2 ) 

dx 2 -— • 

Xi - X 2 + It 


(37) 


In order to obtain the dT /(/.'^-contribution, we have to im¬ 
pose the representation (19 1 on B^ 2 \x i,x 2 ) and, then per¬ 
form the integration over dx 2 by part. However, as shown 
above, B^ 2 \x 1 , x 2 ) does not have the representation (19 1 . 


This property seems to be natural from the point of view of 
gluon poles relation m to Sivers functions as the latter is re¬ 
lated to the projection 7 + . As for the function _B^(x 1 , x 2 ), 
the transverse derivative of Sivers function resulting from tak¬ 
ing its moments may act on both integrand and boundary 
value. Our result suggests that only the action on the bound¬ 
ary value related to B^\x i,x 2 ) should produce SSA. It is 
certainly not unnatural keeping in mind that the integrand dif¬ 
ferentiation is present even for simple straight-line contours 
which are not producing SSA. 
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FIG. 1: The Feynman diagrams which contribute to the polarized Drell-Yan hadron tensor. 




FIG. 2: The matrix element (correlator) of nonlocal twist-3 quark-gluon operator within the momentum representation. Here t = ki — k\ and 

fci = k 2 = {x 2 pt, , k 2 x) 
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